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Atherosclerosis is a progressive disease characterized by the thickening, hardening and 
loss of elasticity of inner arterial walls. The pathologic process of this disease underlies most 
coronary heart disease (CHD) and strokes. Atherosclerosis has plagued Western societies and 
is responsible for over half of the deaths among the elderly worldwide. As a result, the disease 
has been the subject of research for over two centuries [1-37]. Research has intensified 
worldwide because atherosclerosis is becoming the biggest killer among otherwise healthy 
individuals. The World Health Report estimated that atherosclerosis would be the greatest 
single disease burden worldwide by 2020 [38]. The direct and indirect costs associated with 
the disease have been estimated at SI 1 1.8 billion every year in the United States alone [39]. 
Despite our growing knowledge of the cause of atherosclerosis, the prevalence of the disease 
has yet to be reversed. 

Over the past five decades, lipid-lowering therapy has dominated strategies for 
preventing and treating atherosclerosis [20,40-43,68,69]. The lipid hypothesis suggests that 
the elevated level of low-density lipoprotein (LDL) or cholesterol in human blood is the 
primary cause of the disease. It is reported that there are about 55 million American adults 
who have elevated levels of LDL that warrant intervention [44]. However, clinical evidences 
indicate that many individuals in the United States develop CHD in the absence of 
abnormalities in lipoprotein profile [37]. Recently, a number of studies focused on 
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inflammation as a primary cause in atherosclerosis [37,45-48]. The inflammation hypothesis 
emphasizes that macrophages have essential functions in all phases of atherosclerosis. 

From epidemiological and clinical studies, numerous factors have been linked to the 
disease. Contributory factors that have been linked to atherosclerosis include elevated LDL 
level, reduced high-density lipoprotein (HDL) level, hypertension, smoking, family history, 
systemic inflammation, infectious agents, age, high-fat diet, nutritional status, immunity, 
obesity, hormones, lack of exercise and emotional stress [20,32]. Additionally, clinical and 
experimental evidences have associated elevated heart rate with the disease [49,50]. However, 
the presence of any combination of these factors will not guarantee atherosclerosis nor will 
their absence necessarily indicate prevention from the disease. Furthermore, these factors 
affect the entire body and cannot account for the extreme localization associated with 
atherosclerotic lesions. 

The objective of this chapter is to resolve some of the fundamental challenges in studying 
atherosclerosis by creating the Mass Transfer Flux Model and the Lesion Adherence Model 
using analytical methods. The Mass Transfer Flux Model is comprised of a bioheterogeneous 
reaction model, a natural convection model and a boundary value model, and is used to unite 
atherosclerotic risk factors with potential clinical applications. The Lesion Adherence Model 
is comprised of a thermodynamics model and a boundary value model, and is used to further 
the understanding of atherosclerotic mechanisms. We aim to show that analytical methods 
can play a crucial role in the study of atherosclerosis when the reward of prevention and the 
potential for regression is maximal. 

1. Modeling Atherosclerosis 

As with many major human diseases, the mechanisms involved in atherosclerosis are 
complex and multifaceted. The characteristic evolution time of this progressive disease is 
slow, and clinical symptoms only manifest after years of development. These conditions 
make atherosclerosis especially difficult to study in vivo. For these reasons, we must rely on 
models of the disease's processes. Atherosclerotic models permit the simulation of cause and 
effect which is used in place of a real atherogenic event to draw inferences regarding the 
response of the event to changes in various input conditions. 

Take an atherogenic process to be the first diagramed box in Figure la. Inputs to this box 
are determined by environmental features including interactions among the blood constituents 
such as lipoproteins, monocytes, and platelets, characterization of local hemodynamic flow, 
and reactionary responses by arterial endothelium cells at the plasma-endothelial interface. 
The output corresponds to the various manifestations of the disease that result from the 
atherogenic process with the given the input conditions. These manifestations may include 
lesion growth rate, size, shape, composition, location, and proclivity to rupture. The 
manifestations may be regarded as a collection of ultimate clinical goals that include 
prediction of disease risk level and primary cause of disease, establishment of therapeutic 
treatment, and assessment of treatment efficiency. 
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Figure 1 : The philosophy of modeling atherosclerosis 

To reveal the atherogenic event, we replace atherosclerosis with a physical system 
(Figure lb) in which the system is defined for some definite region of space. We then assume 
that the outputs from the system are approximately the same as the outputs that would have 
resulted had the real event been in the first atherosclerosis box. All the features of the real 
atherogenic event will not be present in the physical system, but using a reliable hypothesis 
for the atherogenic event and a set of efficient parameters as inputs for the system, we can 
equate the two box processes. 

In order to truly model a process, we need to replace the physical system with a series of 
mathematical models (Figure 1c) by imposing entity balance and conservation principles of 
mass, momentum and energy to the physical system and creating a set of mathematical 
variables as inputs to the model. The outputs of the mathematical model are the ultimate goal 
of modeling and are determined by solving a system of differential equations. We will take 
the above-mentioned approach to modeling atherosclerosis. 



Entity Balance 

The entity balance of a general system is expressed as 

Input + Generation = Output + Accumulation ( 1 ) 

in which the input is defined as quantities cross the system boundary from the outside to the 
inside in a time interval A/, and the output is defined as the converse. The accumulation is 
defined as the difference between the starting quantities and the ending quantities within the 
system in this time interval. The generation is defined as the quantities that are produced 
within the system in this interval. Conservation principles of mass, momentum and energy 
govern the transport processes without either creating or destroying the quantities. The entity 
balance has the form 



Input = Output + Accumulation 



(2) 
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In a steady-state process, the quantities do not change with time. For processes of rates, the 
entity balance has the form 

Input Rate + Generation Rate = Output Rate (3) 

Continuity assumes that the fluid has a unique value at a given point and time. We will 
consider our models to be a conservative, steady and continual system. 

2. Theories of Atherosclerosis 

Many theories for the origin and development of atherosclerosis have been postulated 
over the years. However, none of them gives a complete and satisfactory picture of the 
mechanism of the disease. Several major theories are summarized below, but current views 
are far less categorical and usually incorporate elements from any number of these 
hypotheses. 

Thrombogensis 

This theory states that mural thrombosis is the essential and primary step in atherogenesis 
and was first proposed by Rokitansky in 1 852 [2]. After thrombi forms in the vessels, they are 
incorporated within the intima. A sequence of reactions including the migration and 
proliferation of smooth muscle cells takes place, which finally leads to the atherosclerotic 
lesion [12]. 

Lipid Theory 

This theory regards atherosclerosis as a consequence of defects in the metabolism of 
cholesterol and other lipids and was first noted by Anitschkow in 1913 [4]. The theory points 
to high serum cholesterol or LDL levels as a major risk factor of the disease [20,40-43]. Low- 
density lipoprotein (LDL), a major cholesterol-carrying lipoprotein in the serum, has been 
recognized by the National Cholesterol Education Program (NCEP) as the major atherogenic 
lipoprotein [43]. 

Inflammation Theory 

This theory was first proposed by Virchow in 1856 [3] and emphasizes an early event in 
atherosclerosis as an inflammatory response to an injured arterial wall. Based on the 
inflammation theory, Ross and Glomset in 1973 proposed the response-to-injury hypothesis 
[10]. The hypothesis is based on the idea that the functioning of the endothelium is changed 
locally due to mechanical damage imparted on the endothelium. These changes lead to 
secretion of platelet-derived growth factor, migration and proliferation of smooth muscle cells 
in the intima and further endothelial damage. As a consequence, atherosclerosis arises as an 
exaggerated response to the initial injury to the endothelial layer of the vessel wall, which 
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leads to an atherosclerotic lesion. Current research emphasizes macrophages as ha\ 
essential roles in all phases of the disease [37,45-48]. 



\ 



Flow Theory 

This theory relates atherosclerosis to the effects of the bloodstream on the arterial wall. 
Experimental findings reveal that lesions are often found near branch sites, curved vessels, 
bifurcations of large arteries, and other regions of perturbed blood flow. A number of 
hypotheses arise from the flow theory [6-8, 1 3- 14, 1 8,2 1 ,35,5 1 -56] : 



Stagnation Point Hypothesis 

Proposed by Fox and Hugh in 1966 [6], this hypothesis states that the localization of 
atheroma is characterized by boundary layer separation where secondary flow and blood flow 
perpendicular to the wall may occur. According to this theory, stagnation points where blood 
components spend longer time than normal or where the flow causes increased transport of 
blood components may be the cause of atherosclerosis. 

High Wall Shear Stress Hypothesis 

Proposed by Fry in 1968 [7], this hypothesis focuses on the notion that shear forces 
caused by blood flow passing by the endothelial wall may be so high that the vessels become 
damaged and the subendothelium exposed, thus leading to an atherosclerotic event. 

Low Wall Shear Stress Hypothesis 

Caro et al in 1969 [8] proposed that atherosclerotic lesions occur in areas experiencing 
low and fluctuating wall shear stress. This hypothesis views mass diffusion as the mechanism 
of LDL transport to the vessel wall, which increased at low local shear rates. 

Diminished Lateral Pressure Hypothesis 

Texon in 1980 [13] developed a concept called the hemodynamic basis of atherosclerosis. 
This hypothesis states that atherosclerotic lesions occur at the segmental zones of diminished 
lateral pressure generated by the flowing blood in the region of arterial bifurcations, 
curvatures or branching called the lesion-prone sites. 



Convection-Diffusion Hypothesis 



When endothelial wall lining is damaged, platelets adhere to the subendothelium and 
form thrombi. This theory focuses mainly on the platelet-vessel wall interaction caused first 
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by convective and diffusive transport of platelets to the vessel wall, then adhesion to the 
vessel, and finally aggregation of platelets to form a thrombus [71]. 

In devising these theories of atherosclerosis, three well-established investigative methods 
have been used. The oldest method relies on many valuable observations of atherosclerotic 
plaques in autopsies within various age groups and with diseases such as hyperlipidemia, 
diabetes and hypertension. The second method is the use of epidemiological studies to explain 
factors that promote or protect people from atherosclerosis. The third method utilizes classical 
experimental pathology to study sequences of lesion development or regression in animal 
models of atherosclerosis. 

While each of the above theories has its own substantiating clinical or experimental • 
evidences, they are not mutually exclusive, and they only attempt to explain a part of the 
many aspects of atherosclerosis. Even though these theories and methods are worthy of 
investigation and each one has contributed to an increased understanding of atherosclerosis, 
little substantial progress has been made in the prevention or regression of the disease [38]. 
Atherosclerosis is a multifactor disease with many risk factors. Analytical models are useful 
in studying atherosclerosis because the models can explore interactions among numerous 
factors at the same time, thus providing a better description of the disease. 

3. Biological Environment of Atherosclerosis 

The hypothesis of atherogenesis relies heavily on the understanding of atherosclerotic 
mechanisms and the evolution of atherosclerotic lesions. Atherosclerotic events may be 
considered to be biological responses of endothelial cells to blood constituents and local 
blood flow on a plasma-endothelial interface at the arterial wall. Thus, we must first consider 
the biological factors and their manifestations in atherosclerosis. 

3.1. Blood Constituents 

The blood is a complex suspension of red and white corpuscles, platelets and lipid 
globules dispersed in a colloidal solution of proteins. The blood is roughly 45% cells, the 
largest of which are red blood cells (-5 x 10 6 per mm 3 ). The non-cellular fluid is the plasma, 
which consists of 95% water and 5% proteins. Main players of atherosclerosis in the blood 
include the following: 

Low-Density Lipoproteins (LDL) 

LDL is a major carrier of cholesterol in the blood. High serum LDL has been identified to 
be a major risk factor of atherosclerosis [20, 40-43,68,69]. The sizes of LDL particles range 
from 21 to 25 nm in diameter [57]. LDL size is significant to atherogenesis because LDL 
particles penetrate through endothelial cells and accumulate in the subendothelium where 
they can initiate an atherosclerotic event [20]. In contrast to LDL however, high-density 
lipoproteins (HDL) seem to provide protection from the disease. 
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Monocytes \ 
The circulating leukocyte-monocyte is the first cell to adhere to the intact endothelium. 1\ 
eventually migrates between endothelial cells, locates itself in the subendothelial space, and\ 
transforms into the foamy macrophage through the ingestion of lipids [37,47]. The binding 
and recruitment of circulating leukocyte-monocytes to the arterial endothelium is considered 
a fundamental step in developing an atherosclerotic lesion, in particular those that go on to 
develop plaque rupture. Ruptured lesions typically have a large necrotic core and a disrupted 
fibrous cap infiltrated by macrophages. A number of research studies [37,45-48] suggest that 
monocyte-macrophages have essential functions in all phases of atherosclerosis, from 
development of fatty streaks to processes that ultimately contribute to plaque rupture and 
myocardial infarction. 

Platelets 

The minute, nonnucleated, disk-like cytoplasmic body functions to promote blood 
clotting. Platelets are the source of numerous pro-inflammatory mediators and can initiate the 
vascular phase of an acute inflammatory response. 



3.2. Arterial Wall 

The arterial wall is composed of a number of successive concentric layers. The tunica, 
intima is the innermost arterial layer. It contains a monolayer of endothelial cells lying on a 
basal lamina. The endothelial cells are aligned according to the streamlines in the blood 
vessel and are in constant contact with the blood. At areas of irregular blood flow of lesion- 
prone sites, e.g. bifurcations, the endothelial cell shape is irregular. High shear regions tend 
have more endothelial cells that are elongated whereas low shear regions tend to have rounder 
endothelial cells. Most of the resistance to blood flow in the artery is due to the endothelial 
layer of the intima. A thin sub-endothelial layer, which includes longitudinal collagen fibers, 
elastic fibers, and smooth muscle cells, follows the basal lamina. The internal elastic lamina 
separates the tunica intima with the tunica media. The tunica media is composed of 10 to 60 
layers of smooth muscle cells surrounded by a fine network of collagen and elastic fibers. The 
external elastic lamina separates the tunica media and the tunica adventitia. The tunica 
adventitia is composed of elastic fibers that intermingle with compact collagen fibers. The 
inner elastic fibers running longitudinally while the outer elastic fibers run circumferentially. 
The outermost region of the adventitia is composed of loose connective tissue. 



3.3. Lesions Types 

To successfully model atherosclerosis, we must understand the morphology of the 
disease. Atherosclerotic lesions are never constant but are subject to change as the disease 
progresses through different stages. The following stages defined by the American Heart 
Association (AHA) are used to categorize the evolution of atherosclerotic lesions [1]. 

Isolated Foam Cellsh 

The earliest stage is the intimal thickening due to macrophages that are overloaded with 
lipid droplet inclusions called foam cells. Groups of foam cells are more frequent in eccentric 
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thickening while isolated foam cells are more common in diffuse thickening. The foam cells 
occur in half of infants in the first seven months of life. 

Fatty Streaks 

Fatty streaks start from simple deposition of low-density lipoproteins (LDL) onto the 
aortic endothelium. LDL is transported to the intima through the endothelium. The LDL 
undergoes oxidative modifications to become Ox-LDL, which is potentially cytotoxic. 
Meanwhile, monocytes from the blood stream migrate to the intima and become activated 
macrophages. These macrophages, as well as migrating smooth muscle cells from the media, 
take up the Ox-LDL in large amounts and become the lipid-laden foam cells of the first stage. 
The conglomeration of foam cells, smooth muscle cells, and extracellular tissue is the fatty 
streak and has a yellow appearance due to the accumulation of lipids. Fatty streaks occur at a 
universal age of ten, but progress to more advanced stages at different rates in different 
people. 

Preatheroma 

Preatheroma is a link between the fatty streak and the atheroma and occurs in the middle 
of the second decade of life and only in eccentric thickening. The preatheroma is 
characterized by multiple small pools of particles clustered in the deep intima layer that 
separate and displace some smooth muscle cells. 

Atheroma 

An atheroma is an eccentric thickening that metamorphosed into an advanced • 
atherosclerotic lesion due to increased accumulation and retention of lipid. The extracellular 
lipid is concentrated at the core and replaces much of the musculoelastic intima layer. The 
macrophage foam cells and lipid-laden smooth muscle cells are layered above the core. 

Fibrous Plaques 

Fibrous plaque, also called atheromatous, fibrofatty or fibrolipid plaque, is the 
characteristic plaque of more advanced atherogenesis. The plaque consists of a mass of lipid 
encapsulated in the intima, separated from the lumen by a cap of fibrous tissue. This lesion 
contains not only more smooth muscle cells, but also more collagens, elastic tissues, 
mononuclear phagocytes, and intracellular and extracellular lipids. This is also the stage of 
the disease characterized as arteriosclerosis or arterial hardening because much of the 
elasticity of the inner arterial walls is lost. 

Complex Lesions 

The final stage is the complex lesions that are of thrombus formation with deposits of 
fibrins and platelets. Complex lesions no longer have the endothelium covering and blood 
platelets are able to adhere to the subendothelial tissue leading to the formation of a thrombus. 
Complex lesions may hemorrhage and have the same white appearance as the fibrous plaques. 

3.4. Hypothesis of Atherosclerosis 

From the previous discussions, we propose that the mass transfer flux of LDL and 
monocytes in the blood to the arterial endothelium at lesion-prone sites is a primary cause of 
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atherosclerosis. We will motivate this hypothesis by using analytic methods presented in th 



The second phase of modeling atherosclerosis is to translate the biological environment 
and parameters that affect the disease into physical parameters founded on physical principles 
like fluid dynamics and thermodynamics. In particular, we aim to model two processes in 
atherosclerosis: (1) the transport of atherogenic constituents such as LDL from the blood to 
the arterial endothelium and (2) the interaction of transferred constituents with the. arterial 
wall which causes lesions to form A bioheterogeneous model and a natural convection model 
describing mass transfer flux of atherogenic constituents will be used to model the first 
process. A thermodynamics model describing lesion adhesion will be used to model the 
second process. Some scientific terms and quantities are employed for the physical system, 
and they must first be properly defined and understood. Refer to references [58-62] for further 
information of these terms and quantities. 

4.1. Hemodynamic Behavior 

Under certain physiological conditions, blood exhibits laminar flow in which parallel 
streamlines of fluid governed by viscous forces slip by one another in an orderly manner. In a 
circular tube such as the artery, laminar flow may exhibit a parabolic-shaped profile. The 
center of the laminar flow profile has the maximum velocity. At the wall of the tube, the 
velocity of the flow is zero meaning that the layer of fluid next to the wall remains fixed. This 
condition of no velocity at the wall is also referred to as the no-slip condition. 

Under other conditions, blood flow changes from an orderly streamlined laminar flow to 
a disorderly turbulent flow. This occurs when the velocity field of the blood fluid is 
stochastic. Compared to laminar velocity profile, the turbulent flow velocity profile is much 
blunter at the central portion of the tube and much sharper at the walls. In addition, turbulent 
flow tends to form a number of eddies. Atherosclerotic lesions are often found at sites of 
turbulence in the arterial system because the mass transfer flux of atherogenic particles rely 
on turbulence eddies. 

The criterion to distinguish a flow as turbulent or laminar is the dimensionless Reynolds 
number. The Reynolds number, Re/\s d-u Q /v where d is the diameter of the blood vessel, 
u 0 is the mean blood velocity, and v is the kinematic viscosity of the fluid. Re is physically 
interpreted as the ratio between inertial and viscous forces. At low Re, blood flow is 
dominated by viscous forces which result in orderly laminar flow. When inertial forces 
dominate the flow, Re is high and the flow is turbulent. In general, blood flow is laminar at Re 
less than 2000. The critical Re is the transitional limit of the flow from laminar to turbulent. 
For transition to turbulence, blood flow through an artery generally has a critical Re of -2300. 
In a bifurcation, however, critical Re is -600 and can be as low as 400 for flows with high 
risk of turbulence [55]. 
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The boundary layer is a crucial element of blood flow in the arteries because exchange of 
cells or blood constituents, e.g. LDL, between the blood and the arterial wall is dominated by 
convective and diffusive phenomena within the boundary layer. Boundary layer theory is 
based on the assumption that close to the vessel wall there exists a thin layer of fluid that is 
mostly governed by viscous forces and outside that layer, those viscous forces may be 
neglected. In the boundary layer, the velocity changes from zero at the wall to a velocity 
unimpeded . by viscous forces. The boundary layer remains intact unless the angle of fluid 
motion along the wall surface is too large in which case the boundary layer may detach itself, 
creating vortices and wakes. Plasma flow in the boundary layer near the arterial wall is 
laminar. 

Finally, viscous properties of blood can also affect its hemodynamic behavior. The 
viscosity of fluid is defined by /j-t I y where // is the viscosity; r is the shear stress and y 
is the shear rate. A fluid is known as a Newtonian fluid when its viscosity at a specified 
temperature is constant. Plasma behaves generally as a Newtonian fluid. The whole blood, 
however, is non-Newtonian because of its aggregation properties at low shear rates, but it 
may be considered reasonably Newtonian for shear rates greater than 100 sec x and vessel 
diameter greater than 0.3 mm. 

4.2. Heterogeneous Reaction 

Heterogeneous reaction is an important physicochemical operation. The reaction' takes 
place at certain phase boundaries and is greatly determined by hydrodynamic factors. This 
reaction consists of three stages. The first stage is the mass transfer flux of reaction particles 
in a solution to a reaction surface at the phase boundaries. The second stage is the 
heterogeneous reaction itself at the surface. The final stage is the reaction products that depart 
from the reaction surface. Our models will focus on the first process, the mass transfer flux 
stage, because both second and third stages are dependent on the first. 

4.3. Mass Transfer Flux 

The transport of blood solutes such as LDL from the blood to the arterial wall is a major 
contribution to atherosclerosis. The mass transport of these blood solutes is governed by 
several different mechanisms. The first of these mechanisms is the biomacromolecular 
diffusion of solutes induced by a concentration gradient. Consider a low solute concentration 
in human blood with a constant diffusion coefficient. The diffusive flux has the form 

J c =-DVc (4) 

where D is the diffusion coefficient, Vc is the gradient of the solute concentration 

distribution, and J c is the diffusive flux in the direction of decreasing concentration. The 

second mechanism is the convection of the solute in which circulating LDL are transferred 
through the bloodstream. The convective flux has the form 
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where u { is the velocity vector of blood fluid. The third mechanism, natural convection, is 
caused only by the density gradient of the solute in a fluid layer near the reaction surface. The 
gradient is derived from the concentration change of LDL at the surface. The mass transfer 
flux of blood solutes, such as LDL, from the bloodstream to the arterial endothelium of 
lesion-prone sites is governed by these three mechanisms. 

4.4. Thermodynamics Interface 

Atherogenesis is related to nanoscale fluid dynamics and macromolecular transport at the 
arterial endothelium [29]. The behavior of atherosclerosis is time-dependent, and the 
dominant phenomena may occur in the range of 0-100 A. Thus, infonnation about the surface 
is very difficult to obtain. The adherence of an initial lesion on the arterial endothelium can be 
treated as a thermodynamics problem in a mixing system because thermodynamics deals with 
macroscopic and statistical behaviors of interface rather than with details of their molecular 
structures. The terms surface free energy, work of adhesion, surface stress, and contact angle 
are used in describing a plasma-endothelial interface in the thermodynamics model. 

The surface free energy, y , a major thermodynamics quantity used to characterize an 
interface, is the reversible work used to create a unit area of surface at constant temperature, 
volume, and chemical potential /(//,). The surface free energy y of a newly created surface is 
the amount of surface work (dw) done to form the new surface area (dA) and is defined as 

r = ^ (6) 
' dA 

In liquid, but not in solid, the surface free energy can be interpreted as the surface stress or 
tension, which is the work necessary to stretch or compress an existing surface. 

When phase A adsorbs to phase B, there will be a surface energy of adhesion, E AB , or 
also called the work of adhesion, W AB , between the two different phases. The Helmholtz 
interfacial free energy (or tension), y AB , is defined as 

A newly formed surface may restructure with time and will reach some equilibrium with its 
surrounding environment. 

According to Young's equation, the equilibrium of the interfacial tensions at a liquid- 
solid interface is expressed as 



(8) 
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where 9 is the equilibrium contact angle, y s is the solid surface free energy, y L is the liquid 

surface free energy and y SL is the surface free energy at the liquid-solid interface. Combining 

Eq. (7) and (8) results in the work of adhesion, W SL , which is expressed using the contact 
angle 

W sl =y l {\ + qos0) (9) 

The contact angle technique is one of the most sensitive methods known for obtaining true 
surface information. 



5. Mathematical Models 

Since the blood flow is confined to a circulation system, the flow obeys the conservation 
principles of mass, momentum and energy. The entity balance and conservation principles are 
used to construct the mathematical models. 



5.1. Continuity Equation 

Consider the mass M of a blood density p(x,t) in a volume V 

M= [p(x,t) dV (10) 



where p{x,t) is a continuously differentiable function of spatial coordinates x(x 1 ,x 2 ,x 3 ) 
and time /. Given that masses are neither created nor destroyed in a conservative system from 
Eq. (2), we impose the condition DM I Dt = 0 to yield the continuity equation of the blood 
fluid 

^ + V.(pi/ # ) = 0 (11) 

where u^xj) is the velocity vector with components (u l ,u 2J u 3 ) and V is the gradient 
vector operator. 



5.2. Equations of Motion 

The conservation of momentum states that the rate of change of momentum of the body 
is equal to the force acting on the body, which is expressed by the following equation of 
motion 
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at 



(12) 



where cr.. is the stress tensor where ij = 1, 2, 3 and F,(x„/) denotes the body force 
vector per unit mass. The stress tensor of the fluid is (J.. = -p,I ¥ + d.. where - is the 
isotropic stress tensor with unit tensor I i} at a hydrostatic pressure p, and d.. is the non- 
isotropic stress tensor or the deviatoric stress tensor caused entirely by the fluid motion. The 
fluid is assumed to be isotropic and d y is assumed to be linearly dependent on the local 
velocity gradients. The fluid under these two assumptions is known as a Newtonian fluid with 
the following constitutive relation 



d, = 2ft 



U = 1,2,3 



(13) 



where e ti - — 
9 2 



dxj dx i 



is the symmetric part of the velocity gradients tensor called the 



rate of strain tensor, and fi is the fluid viscosity. 

Substitution of Eq. (13) into (12) yields the Navier-Stokes equation that governs the 
motion of fluid 



= v 



2n[e ij -\{V-u i )l 9 ] 



(14) 



If we take the temperature to be uniform over the fluid, an isothermal fluid, then the Navier- 
Stokes equation reduces to 



— '- + (u r V)u, 
ot 



v 2 „.+lv(v-«.) 



-Vpt+pFt 



(15) 



where V 2 denotes the Laplace operator. For an incompressible, isothermal and Newtonian 
fluid, Eq ( 1 1 ) and ( 1 5) further reduces to 



(16) 



dt 



= pF i -V Pi +//V 2 M; 



(17) 
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du. 

where p — L is the transient acceleration term, p(u i • V)w. is the convective acceleration 
dt 

term, pF i is the body force term, — Vp. is the pressure gradient term, and juV 2 u i is the 
viscous force term. 

For an axi symmetric flow in a tube, Eqs. (16) and (17) is expressed in cylindrical 
coordinates as 



du r u r du 7 
— ^-h-^ + — ^ = 0 
dr r dz 



(18) 



du z du z du, \ dp fi 

— - + u r — - + u z — - = F z + — 

dt dr dz p dz p 



( d 2 u z 1 du z d 2 u 

7- + L + ; 

^ dr 2 r dr dz 2 



2.. ^ 

z 



(19) 



du 



du. 



^- + u r — L + « z — = F r - + 

dt r dr z dz r pdr p 



1 dp p( d 2 u r 1 du r u r d 2 u 



dr 2 r dr 



dz 2 



(20) 



where u r is the radial velocity component and u. is the axial velocity component. 



5.3 Convective Diffusion Equation 

Consider a conservative system of volume V that contains blood particles. A surface S 
bounds the volume. Following mass conservation and entity balance, the transport of particle 
mass without creating or destroying the quantities in the system can be expressed as 



ht j 



(21) 



where J = cu i - DVc is the sum of diffusive and convective flux from Eq. (4) and (5). 

The left side of Eq. (21) represents the rate of change of the number of particles in the volume 
and the right side is the number of particle that moves into the volume. Eq. (21) can be 
rewritten in the differential form as 

dc 

— = -div (D grad c) - div (c u i ) (22) 
dt 



For an incompressible fluid with a constant diffusion coefficient D, Eq. (22) simplifies to the 
following .convective diffusion equation 
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— + («, V)c = £)V 2 c (23) 
ot 

Eq. (23) in Cartesian coordinates is 

dc dc dc dc ^,d 2 c d 2 c 3 2 c v 

— + w — + v — + w — = £> ( — - + — - + — -) (24) 
dt dx 3y dz V dy 2 dz 1 ' K } 

6. Analytical Models 

6.1. Mass Transfer Flux (MTF) Model of Atherogenic Particles 

Following the terms, concepts, and methods in the Physical Systems and Mathematical 
Models sections, we create a bioheterogeneous reaction model, a natural convection model 
and a boundary value model to investigate the mass transport flux of LDL and monocytes 
from the blood phase to the atherosclerotic lesion phase. The symbols used are denoted in the 
nomenclature at the end of the chapter. 

Bioheterogeneous Reaction Model 

Bioheterogeneous Reaction 

The bioheterogeneous reaction of atherosclerosis occurs at a plasma-endothelial interface 
and is largely dependent on local hemodynamics. This reaction consists of three stages: The 
first stage is the mass transfer flux of LDL and monocytes in blood to a reaction surface of the 
arterial endothelium at the lesion-prone sites. The second stage is the bioheterogenous 
reactions between these transferred solutes and endothelial cells at the surface, which leads to 
LDL particles and monocyte cells being filtered onto the surface and accumulating in the 
subendothelium. The final stage is the promotion of atherogenesis by the reaction products at 
these sites. We will focus on the first mass transfer flux stage because both the second and 
third stages are dependent on the first. 

Local Blood Flow 

Clinical and experimental observations [63-67] show that atherogenesis is often found in 
regions of arterial bifurcations, branching or curvatures, where the bloodstream is 
characterized by disturbed non-axial flow, reverse flow, flow separation and stagnation. 
Following Sections 4.1 Hemodynamic Behaviors and 4.3 Mass Transfer Flux, the blood flow 
in regions of lesion-prone sites may be nominally divided into three zones. The first zone is 
the disturbed flow region farthest from the reaction surface. In this zone, both solute 
concentration and density are constant because blood constituents and momentum are 
transferred by eddies. The second zone, nearer to the surface, is called the momentum 
boundary layer where the plasma laminar flow by molecular viscosity occurs. In the boundary 
layer, there are solute concentration changes and density gradients caused by the solute 
concentration change at the surface. Since kinematic viscosity of fluid is larger than its 
diffusion coefficient, the momentum and solute in the layer are mainly transferred by eddies 
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that rapidly decay toward the surface. Finally, the third zone, called the diffusion boundary 
layer, is nearest to the reaction surface where eddies are negligible, and LDL or monocytes 
are transferred by diffusion. The fluid motion in the layer is governed by natural convection. 

Consider blood flow through a branched vertical artery tube of internal radius a,, where 
(r, 0, z) is a set of cylindrical coordinates with the z-axis coinciding with the axis of the 

tube. The flow at the border between the momentum boundary layer and the disturbed flow 
zone has to satisfy the following continuity condition: 

u r =V d at r = a-h # (25) 

where u r is the radial velocity component of plasma fluid in the momentum boundary layer, 

V d is the eddy in the disturbed region and h is the thickness of the momentum boundary 

layer. Effects of local hemodynamic factors on the bioheterogeneous reaction are determined 
by the interactions of blood flow among the three mentioned zones. 

Natural Convection Model 

Governing Equation 

Using same notations, we investigate the plasma flow over a reaction surface in the 
momentum boundary layer in the range of {a - h) < r < a , 0 <0 <2n , - oo < z < oo 
and under body forces because atherosclerotic lesion arises from the layer. Denote the 
concentration of solute far from the surface as c 0 and at the surface as c, . As indicated in the 

second stage of the bioheterogenous reaction (Section 4.2), a decrease in c x corresponds to 
LDL or monocytes being filtered into the subendothelium. If the concentration difference 
( c 0 — C, ) is sufficiently high, natural convection in a fluid layer near the surface will be 

generated with a velocity component u z in the axial direction of the tube, which will vary 

with the radius r . The natural convection model corresponds to a physical situation of 
atherosclerotic event at the location. Since plasma fluid is known to be Newtonian, we will 
use the set of differential equations (11), (17) and (23) outlined in the Mathematical Models 
section to investigate the physical situation at the momentum boundary layer, 
Assumptions 

Based on the Physical System section, the following assumptions are used to simplify 
Eqs (11), (17), and (23): 

(1) The plasma is an isothermal and incompressible fluid. 

(2) The rate of plasma flow is steady. 

(3) The pressure gradient in the momentum boundary layer is negligible because the 
pressure field in the disturbed flow zone is constant. 

(4) The solute concentration is low and the rate of concentration change is steady. 

(5) The diffusion coefficient is considered constant due to low solute concentration. 

(6) The viscous dissipation in the fluid is negligible. 
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(7) Edge effects become negligible since the length of the tube is long compared to the 
region being studied. 

(8) The flow does not slip on the endothelial surface of the arterial vessel. 

(9) The diameter of the tube does not vary with internal pressure. 

These assumptions are used in most studies of hemodynamics and mass transport 
[13,54,58,60,61]. 

Simplified Governing Equation 

Because the momentum boundary layer overlaps the diffusion boundary layer in the 
model, we analyze the plasma flow under body forces in the momentum boundary layer. 
Since a » h , the flow is considered planar in a vertical semi-infinite plane at 0 < y < h 
and - oo < x < oo , where x-axis is along the plane surface and toward the flow direction, and 
jK-axis is outward normal of the surface. Using Boussinesq's approximation that neglected 
density variations with concentration in all terms except for the body force, Eqs. (11), (17), 
and (23) in Cartesian coordinates are simplified to: 



du y du v 



+ 



dy dx 



(26) 



fox du x Fx d 2 w r 

u v — - + u r — - = — + v ^ (27) 

y dy x dx p dy 2 1 1 

dc dc ^ d 2 c 

''H*'-Sl- D W (28) 

where u x and u y are the velocity components of fluid in the x and y directions respectively, 
F x is the body force component consisting of the transient inertial force and gravitation in 
the x direction and v = fit p is the kinematic viscosity of plasma. Eqs: (27) and (28) have 

u - ,t^w t • * d2u * Qlu x * 9 2 c d 2 c 
been simplified by taking mto account ~ » and » . 

dy 2 dx 2 dy 2 dx 2 

Boundary Conditions 

According to the previous models and assumptions, the boundary conditions of the 
boundary value model are: 

u x =0 > u y =°> c = 0 at 7 = 0 (29) 



u y = v d at y = h 



(30) 
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c = c Q at y -» oo 



(31) 



Mass Transfer Flux 

Applying conditions (29)~(31) to Eqs. (26)-(28) yields the following concentration 
distribution of blood solute in the diffusion boundary layer 



c = 0.69 c 0 y 



(32) 



Substituting Eq. (32) into (4) yields the mass transfer flux of the solute to the reaction 
surface: 



J = 0.69 c n 





27 f 






I J 


V 



(g cos a? + cou^k 



(33) 



with 7 being the flux toward the inner artery wall, g being the gravitational acceleration, 
atherosclerotic parameters c 0 being the LDL or monocyte concentration in human blood, co 

being the heart rate, u 0 being the transient mean velocity of blood fluid in the axial direction 
of the arterial vessel near the studied region, v d being eddies determined by local dynamics 
of the fluid, v being the kinematic viscosity of plasma fluid in the region, x being the 
diffusional length of LDL or monocyte along the inner artery wall, D being the diffusion 
coefficient of atherogenic particles, k being the density parameter determined by 
concentration change and a being the angle between the mean velocity and gravity. We will 
refer to Eq. (33) as the atherosclerotic mass transfer flux equation or AFE for short. 



6.2. Analytical Results and Clinical Implications of MTF Model 

The MTF model provides a means to analytically study the relationship among 
atherogenic factors including blood LDL or monocyte concentration, heart rate, transient 
mean blood velocity, eddy velocity, kinematic viscosity, diffusional length and diffusion 
coefficient of atherogenic particles, the density parameter, and the angle between the mean 
velocity and gravity. These factors are combined through Eq. (33), the atherosclerotic mass 
transfer flux equation (AFE), which describes the total flux of atherogenic particles to the 
lesion prone sites where an atherogenic event can occur. 

Assessing Lipid-Lowering Therapy 

From the AFE, a 1% reduction in LDL concentration results in a 1.2% reduction in the 
disease risk. Epidemiological studies have shown that a 1% reduction in an individual's total 
plasma LDL level lead to a 1.5 to 2.0 % reduction in the risk [68,69]. The clinical studies are 
in good agreement with results of the MTF model and indicate that the model can provide a 
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reliable theoretical baseline for lipid-lowering therapy. By determining the lipid concentration 
change due to therapy, we study the impact of the concentration change on mass transfer flux 
of LDL using AFE and quantitatively assess therapeutic efficacy of the treatment. 
Epidemiological results that show greater reduction in risk than our predicted results may 
indicate additional protection from disease by the therapy. 

Non-independent Contribution of Risk Factors 

While the relative impact of risk factors is determined by the individual's physiological 
parameters, the contribution of risk factors to the disease may vary greatly among individuals 
as predicted by the AFE. For example, autopsy and clinical studies [13-14,16,22] suggested 
that regions of arterial bifurcations, branching or curvatures had the greatest predilection for 
atherosclerosis. Internal angles among 70 human aortic bifurcations can vary widely from 10 
to 70 degrees [70]. Different internal angles in physiologically different bifurcations amongst 
people may lead to different a in the AFE. Thus internal angle may be a greater risk factor 
under some conditions. To illustrate, take two different internal angles 15° and 45°. Under a 
1% increase in individual's LDL level, AFE produces a 7.2% higher risk for 15° than for 45°. 
This 7.2% risk caused by the difference in bifurcation angle is significantly higher than the 
1 .2% increase in risk from 1% increase in LDL level alone as mentioned previously. The case 
shows that arterial geometry in certain instances can play a greater role in atherosclerosis than 
simply LDL level alone. Therefore, AFE not only can give relative impact of various risk 
factors, but can also show how certain risk factors can dominate others depending on the 
person's physiological parameters. Simply put, the contribution of risk factors is conditionally 
dependent on the person's physiology. There is not a simple independent relationship 
between physiological parameters and the risk of disease, but rather the risk of disease is a 
combination of the physiological parameters themselves and their impact on one another. 
This observation may explain the sometimes inconsistent results that are produced from just 
using independent risk factors to explain the cause of atherosclerosis [72]. The MTF model is 
therefore useful in determine the interaction among the atherogenic parameters and their 
overall impact on risk of disease. 

Optimizing Therapeutic Treatment 

Epidemiological studies suggest that elevated LDL level in blood, inflammation, 
infectious agents, smoking cigarette and depression are risk factors [20,31,32]. Both 
experimental and clinical evidences indicate that increased heart rate is also associated with 
atherosclerosis [49,50]. However, it is difficult to assess the combined contribution of these 
factors to the disease. The MTF model can be used to quantitatively assess the risk factors 
and generate an efficient sequence of therapy for the treatment of the disease. 

For example, AFE predicts a 0.2% increase in the risk of disease from a 1% increase in 
heart rate, which may be induced by smoking or stress. A 1% increase in leukocyte-monocyte 
concentration in the blood, which may be cause by infectious agents, results in a 1.2% 
increase in the risk of disease. A 1.0% increase in both the individual's monocyte level and 
heart rate leads to a 1 .4% increase in the total risk. We therefore conclude that the infectious 
agent is our first therapeutic target and smoking is the second target. This example illustrates 
that AFE can be a potential method to optimize therapeutic treatment that ranks the 
contributions of various risk factors to atherosclerosis. 
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6.3. Lesion Adherence Model 

Following the terms, concepts, and methods in the Physical Systems and Mathematical 
Relations sections, we will create a thermodynamics model and a boundary value model to 
explore mechanisms of atherosclerosis, in particular the initial stage of lesion adhesion, and to 
explain clinical results. 

Thermodynamics Model 

Thermodynamics Model of Atherosclerotic Lesion 

The arterial endothelium is always in direct contact with the blood flow and is borne by 
hemodynamic forces. Endothelial cells have synthetic and metabolic capabilities and act as 
selectively permeable barriers for the transportation of macromolecules. The interfacial forces 
acting on the endothelium can change endothelial functions and. structures including increased 
molecular permeability, LDL accumulation, and endothelial cell damage. In this model, an 
initial lesion adhered on the endothelium is investigated using thermodynamics and energy 
conservation principles. 

Let us consider the blood flow through a branched circular arterial tube of internal radius 
a under an axial transient inertial force /, where (r,#,z) is a set of cylindrical polar 
coordinates with the z-axis coinciding with the axis of the tube. Taking a small fluid element 
near the inner tube wall at a branch point where r = a and z = 0 , the kinetic energy of the 

1 2 

element E k is — • p-V -u t , where p is the density of the fluid, V is the volume of the 

element, and u t is the velocity vector of the element. If an initial lesion occupies the element, 
E k becomes the kinetic energy of the lesion. This lesion is characterized by volume V, axial 
adhesive length of the lesion on the endothelium L a , contact area between the lesion and the 
endothelium A, surface free energy y sl , velocity vector w /? and interfacial shear resistance 
F a , which is the force resisting the movement of the lesion (Fig. 2). 

Energy Equations 

The adherence of an initial lesion on a plasma-endothelial interface can be treated as a 
thermodynamics problem in a mixing system because thermodynamics deals with 
macroscopic and statistical behaviors of interfaces rather than with details of their molecular 
structures. Surface energy of adhesion E sJ is the energy that is required to bring LDL from 

the interior of the blood fluid into the arterial endothelium and form an interface. Consider a 
lesion under an axial transient inertial force/(Fig. 2). 
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Figure 2. Axial cross-section of a lesion adhered on the arterial endothelium. 

Based on thermodynamics and energy conservation principles, the energy equations for 
the lesion are created as 



(34) 



(35) 



E k =E 5l :A 



(36) 



where E s/ =y s +y l - y s i according to the Eq. (7) of the Thermodynamics Interface 
section. 

Eq. (34) describes the equivalence between the kinetic energy of a lesion E k and the 
work done by the total interfacial shear resistance F a • A along the adhesive length L a . Eq. 

(35) represents the equivalence between the surface energy of adhesion E sf and the work 
done by F a at a plasma-endothelial interface. Substitution of Eq. (35) into Eq. (34) yields Eq. 

(36) . If E k is greater than E s{ • A , then the lesion will not form because the kinetic energy 
of the lesion E k overcomes the total energy of adherence E sl • A before the adherence of 
the lesion can occur. 

Dynamic Boundary Value Model 



Governing Equations 

Studies of hemodynamics [6,8,13,21,63-67] suggest that atherosclerotic lesions are 
located in regions of reduced shear stress of the fluid, which is often associated with flow 
separation and turbulence. A transient boundary layer, led by turbulent flow, occurs in these ' 
regions. The studies also indicated that blood flow near inner arterial walls consists mainly of 
plasma fluids. 
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In the dynamic model, the plasma viscous flow in the transient boundary layer near the 
inner tube wall at lesion-prone sites is investigated because lesions arise in this layer. The 
plasma flow through a circular tube is governed by Eq. (17). 

Assumptions 

According the terms and concepts in the Physical Systems section, the following 
assumptions are imposed to simplify the boundary value model: 

(1) The gravitational potential in the fluid has a negligible effect because the transient 
boundary layer is thin. 

(2) The rate of flow is steady, and there is a constant pressure flow. 

(3) The flow does not slip at the inner wall of the tube. 

(4) The length of the tube is long compared to the region being studied, so that edge 
effects become negligible. 

(5) The diameter of the tube does not vary with internal pressure. 

(6) The convection in the flow may be negligible in the model. 

These assumptions are valid in most studies of hemodynamics [13,54,58,61]. 

Simplified Governing Equations 

Let us use the same notations and consider an incompressible plasma flow in a circular 
arterial tube with ranges (a - h(z)) <r <a , O<0 <2k , and - oo < z < oo 5 where h(z) 
is the thickness of the transient boundary layer. Following the above assumptions and using 
cylindrical coordinates (r, 0, z) , Eq. (17) is simplified to 



where u z is the axial velocity component of the plasma fluid, r n is the axial viscous shear 
stress of the fluid, / is the axial transient inertial force, and /i is the viscosity of plasma. 

Boundary Conditions 

According to previous assumptions, the boundary conditions of the model are 




(37) 



and 




(38) 



u z = 0 at r = a 



(39) 



trz =F(z) at r = a(40) 
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where F(z) is the shear force acting on the inner wall of the tube. 

If a fatty streak occurs at r = a and z = 0 , we can substitute F{z) = F a into Eq. (40) 
to yield 



where the interfacial shear resistance F a is defined in the previous model. Through Eq. (41), 
we can now join the thermodynamic model to the boundary value model. 

Compatibility Condition 

As stated previously, there is often turbulent flow at the location of atherosclerosis. The 
flow at these sites consists of two different types. The first is the plasma flow in the transient 
boundary layer that is governed mainly by the viscous force. The second is the blood flow in 
the rest of the regions of the tube that is governed by the axial transient inertial force. Hence, 
the viscous force and the transient inertial force must equal at the border of the layer for a 
unique analysis solution. The compatibility condition is created as 



where f h is the axial transient inertial force per unit area. We may take 



Analytical Solutions 

Eq. (37) and Eq. (38) combined with Eq. (39), (40), and (42) yield the following axial 
shear stress of the plasma fluid , axial velocity component u z , and compatibility 
equation: 



T rz =F a zxr = a,z = 0 



(41) 



T rz =fh at'r = a-A(z) 



(42) 



f k =f.{a-k(z)). 




(43) 



2 =x.(^. r2)+ x. (a _ Mz))2 . ta r^ 

4// 2// \r) 



(44) 



(45) 



Substitution of F(z) = F a into Eq. (43) yields 
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♦ 



F a =f- a-h(z) + 



h(z) 2 



(46) 



2a 



Substitution of Eq. (35) into Eq. (46) yields 




h(z) 



(47) 



2a 



J 



We will refer to Eq. (47) as the adhesion energy equation (AEE). 

Effect of Eddies on Atherosclerosis 

Turbulence eddies do not originate in the transient boundary layer, but these eddies may 
enter the layer from the side r < (a - h(z)) . If we consider the effect of the eddies in the 

layer on atherosclerosis, the thickness of the layer h{z) in the region of the lesion-prone sites 

can be expressed approximately as h{z) ~"B-v d (48) 

where eddy velocity v d is a weak function of z and varies with r . Parameter B is 
determined by the axial velocity u b of the plasma fluid at the border of the layer and the axial 
distance z h between the point being studied and the starting point in the layer. 

The transport of LDL and monocytes in blood to a plasma-endothelium interface at the 
lesion-prone sites may be closely related to v d . Eq. (48) indicates that increased h{z) 

involves increased v d , which implies that LDL in the plasma fluid are more likely to be 
drawn into the endothelium. Increased thickness of the transient boundary layer h(z) and 
increased eddy velocity v d in regions of the lesion-prone sites are risk factors in 
atherogenesis. 

6.4. Analytical Results of Lesion Adherence Model 

Using Eq. (47), the adhesion energy equation (AEE), we can obtain representative results 
using some typical values for a human abdominal aortic vessel [58]. Take the mean internal 
radius of the vessel to be a- 0.5 cm , mean transient axial velocity to be w 0 = 15 cm • s~ l , 

mean transient frequency to be a> = 2 s~ l > density to be p - 1 .0 g - cm~ 2 , adhesive length 
to be L o = 0.2 cm , and the range of the thickness of the transient boundary layer 
0 < h(z) < 0.1 cm . Since f = p-o?-u 0 , we get the axial transient inertial force 
/ = 30 g-cm~ 2 s' 2 . 

To analyze the dynamic boundary value model, substitute r = OA cm and the above 
mentioned values into Eq. (45) to produce a relationship between the shear stress of the 
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plasma fluid and the thickness of the transient boundary layer (Fig. 3). For the 
thermodynamics model, substitution of the above mentioned values into the AEE (Eq. 47) 
produces a relationship between the surface energy of adhesion and the transient boundary 
layer (Fig. 4). 
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Figure 3. The relationship between the shear stress of the plasma fluid and the thickness of the transient 
boundary layer. 
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Figure 4. The relationship between the surface energy of adhesion and the thickness of the transient boundary 
layer. 
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8. Conclusions 

We have demonstrated how analytical techniques can be used to investigate atherogenic 
events such as the transfer of LDL particles from the bloodstream to the endothelial wall and 
the thermodynamics mechanism of how the LDL particles adhere to the wall. The following 
conclusions are drawn from the specifics of the two models. 

8.1. Mass Transfer Flux (MTF) Model of Atherogenic Particles 

We propose that the variables in the AFE be considered as atherosclerotic parameters. 
The mass transfer flux determined by the AFE can be used to determine the causes of 
atherosclerosis and their relative contributions and relationship to the risk of disease. The 
AFE has therefore united some crucial factors from main theories of atherosclerosis and 
integrated risk factors [14,20,31,32,37] into a single analytical parameter, the mass transfer 
flux, in order to evaluate their contribution. 

The analytical MTF model can explain some of the following clinical and experimental 
observations. 

(1) Lipid-lowering therapy plays a central role in clinical prevention and treatment of 
atherosclerosis. Our model shows that LDL concentration should be indeed first 
considered clinically because the LDL parameter c o has the largest impact on the 

AFE. We believe that increased LDL and monocytes levels in blood are a primary 
risk factor because it directly leads to an increase in mass transfer flux and thus 
accumulation of the substance within the subendothelium. 

(2) Not any single one, but a combination of the proposed parameters should be 
considered when devising preventive and therapeutic actions against atherosclerosis. 
Relative impact of the parameters to the treatment can be evaluated using the mass 
transfer flux of the AFE. The flux may therefore be used as a theoretical baseline for 
efficacy of therapeutic intervention and possible regression of the disease. 

(3) No endothelial defects seem to occur during formation of early lesions. This is 
because the plasma viscosity is lower at the lesion formation area than at its 
neighboring sites. According to our MTF model, the mass transport of blood solutes 
to the inner artery wall increases with a decrease in plasma viscosity, which also 
corresponds to a lower wall shear stress at atherogenic sites. Thus lesions tend to 
form near regions of low wall shear stress and high mass transfer flux. 

(4) The mass transport of blood constituents exists not only in the diffusion layer but 
also in the momentum boundary layer and disturbed flow region such as at arterial 
branches or bends of lesion-prone sites. 

The conclusions and analytic results of the MTF model are supported by major clinical and 
experimental evidences [8, 1 3- 1 5,20-2 1 ,3 1 -32,36-37,46-50,68-69] . 
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8.2. Lesion Adhesion Model 

We obtain the following conclusions based on the Lesion Adhesion Model and our 
analytical solutions. 

(1) From Eq. (45), Figs. (3) and (4), we find that the shear stress of the plasma fluid r n 
is proportional to the surface energy of adhesion E sl and inversely proportional to 
the thickness of the transient boundary layer h(z) . This means that reduced T n 
results in increased h(z) and reduced E sl . 

(2) Fig. 4 and Eq. (48) indicate that increased h(z) involves reduced E si and increased 
eddy velocity v d . Increased h(z) in regions of lesion-prone sites is a dominant 
factor in the formation of initial lesions because, based on both of the previous 
models, reduced E sl and increased v d lead to atherosclerosis more easily. 

(3) Clinical and experimental results [8,21,65] suggest that atherosclerosis occurs in 
regions of reduced T rz . According to Conclusion 1 and Eq. (48), reduced 
involves reduced E sl and increased v d . This suggests that LDL in the fluid is more 
prone to be drawn to the plasma-endothelial interface and initial lesions form more 
easily in regions of reduced r n . 

(4) Clinical and experimental results [6,8,13,21,63-67] show that locations of 
atherosclerosis are associated with flow separation and turbulence. Based on 
Conclusion 1 and Eq. (48), lesion formation occurs at these locations because it is 
stimulated by increased h{z) and v d . 

(5) Clinical results [8,13,19,21] show that no endothelial defects occur during the 
formation of initial lesions. According to Conclusion J and Eq. (35), reduced 
involves reduced interfacial shear resistance F q , and since endothelial defects are 
closely dependent on F a , lower value of F a in the lesion formation area than those 

of its neighboring sites is therefore not the cause of endothelial: defects initially. 

(6) Clinical results [1,5,13,16,19,20] show that an initial lesion is caused by depositions 
of LDL on the arterial endothelium. According to the thermodynamics model and Eq. 
(36), LDL rather than HDL are drawn from the interior of the plasma fluid into the 
plasma-endothelial interface because the lower density of LDL yield a lower E sl 

than the E sl for HDL deposition. Therefore, reduced E si or less energy required for 
LDL deposition leads to atherogenesis. 

(7) Clinical results [5,13,16,19,29,58] show that initial lesions do not form in veins and 

capillary vessels even though there is a lower r n in veins than there is in arteries. 

This phenomenon can be explained by the lack of a transient boundary layer in veins 
and capillary vessels. 
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The above conclusions suggest that the locations of atherosclerotic lesions are associated 
with regions of increased transient boundary layer h(z) and reduced surface energy of 

adhesion E sl , which lead to atherogenesis in only a small region of the vast systems of blood 
vessels. Analytical results of the thermodynamics and boundary value model suggest that 
increased eddy velocity v d and reduced E sl in these regions of the arterial branch points are 
the main causes for the formation of initial lesions. 

More specifically, increased h(z) in regions of the arterial branch point is responsible 

for increased v d and reduced viscous shear stress . Reduced E sl and increased v d 

contribute to bring LDL from the interior of the plasma fluid to the plasma-endothelial 
interface, which lead to the formation of an initial lesion. LDL then passes through the arterial 
endothelium, accumulates in the subendothelium, and stimulates monocytes to enter this area. 
The monocytes may then engulf cholesterol, and the lesion progresses to a more advanced 
stage. 

9. Mechanism of Atherosclerosis 

Fung [58] suggested that the most intriguing feature of atherosclerosis is its appearance 
only at certain arterial sites such as bifurcations, branching, and curvatures. It is known that 
many factors influence atherogenesis including dieting, environment, elevated level of LDL, 
inflammation such as rheumatoid arthritis, infectious agents such as Chlamydia pneumonia, 
hypertension, cigarette smoking, family history, genetics, depression, diabetes and obesity 
[20,31,32]. However, these factors affect not only these special sites but also the entire body. 
Out of the vast bodies of vessels in the human circulatory system, only a small segment will 
become atherogenic. What is the mechanism behind the deposition of LDL particles and 
monocytes cells onto these arterial walls and not others? The analytic methods discussed in 
this chapter can help explain the mysterious nature of the extreme localization of the disease. 
Atherogenesis appears only at these locations because the mass transfer flux of atherogenic 
particles and eddies in the blood always occur at these sites. There are no lesions in human 
veins and capillaries because the flux and eddies does not exist in theses vessels. Eq. (33) in 
the MTF model shows that the flux is dependent on the eddy parameter. A number of 
experimental studies indicated that the eddies occurred often at these special locations 
[13,22,63-67], confirming the above analytical results. 

10. Remarks 

The analytical methods discussed in this chapter tackles atherosclerosis as a multifactor 
disease with differently combined risk factors dominating at different stages of disease in 
different individuals. The methods approached atherosclerosis as a mass transport process 
from the blood phase to the lesion phase and is supported by clinical and experimental 
evidences [1,5,13,16,19,32,36-37,47]. 

The analytic methods are sound because they are derived from the entity balance and 
conservation laws of mass and momentum, which govern the atherogenic process. The results 
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given by the analytical methods not only impact the understanding of atherosclera 
mechanism but also can be used in clinical practices. These results may allows physician\ 
predict a total risk of the disease, to determine the primary cause of disease, to assess th 
therapeutic efficacy of treatment, . and to optimize the treatment for the individuals who 
require diagnosis, prevention, and treatment of atherosclerosis. Analytical methods may play 
a crucial role when studying environmental influences, dieting, lifestyle, high blood pressure, 
exercise and genetic impacts on atherosclerosis when the reward of prevention and the 
potential for regression is maximal. 



Nomenclature 

a angle between mean velocity and gravity (degrees) 

L axial adhesive length (cm) 

a 

axial distance between point being studied and starting point of transient 
boundary layer (cm) 
u Q axial mean velocity of fluid (cm I s) 

f axial transient inertial force (gl cm 2 -s 2 ) 

f axial transient inertial force per unit area at border of transient boundary 

layer (glcm-s 2 ) 
u axial velocity component of fluid (cm I s) 

axial velocity of plasma fluid at border of transient boundary layer 
(cm I s) 

F i , F e body force vector (g/cm 2 • s 2 ) 

c, c 0 , c ] concentration of solute ( mg /cm 3 ) 

A contact area between lesion and endothelium (cm 2 ) 

P density of fluid (g/cm 3 ) 

D diffusion coefficient of solute (cm 2 / s) . 

X diffusional length of solute (cm) 

v d eddy velocity (cm I s) 

8 gravitational acceleration (cm Is 2 ) 

M heart rate (s~ x ) 

a inner radius of tube (cm) 

F a interfacial shear resistance (g/cm-s 2 ) 

v kinematic viscosity of fluid (cm 2 1 s) 

E k kinetic energy of lesion (g-cm 2 1 s l ) 

J 5 J u , J c mass transfer flux of solute (mg I cm 2 s) 
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p 


pressure vector (g 1 cm • s 2 ) 




radial velocity component of fluid (cm 1 s) 


F(z) 


shear force acting arterial wall (glcm-s 2 ) 


T rz> T 


shear stress of fluid (g/cm-s 2 ) 


E s , 


surface energy of adhesion (gl s 2 ) 


7i 


surface free energy of liquid (gl s 2 ) 


Yd 


surface free energy of liquid-solid interface (g Is 2 ) 


Ts 


surface free energy of solid (gl s 2 ) 


h(z) 


thickness of transient momentum boundary layer (cm) 




velocity components in Cartesian coordinates (cm 1 s) 


u x 


velocity vector (cm 1 s) 




velocity vector of lesion (cm 1 s) 




viscosity of fluid (gl s- cm) 


V 


volume of lesion (cm 3 ) 



References 

[1] Stary HC, et al. A definition of advanced types of atherosclerotic lesions and a 
histological classification of atherosclerosis: A report from the Committee on Vascular 
Lesions of the Council on Atherosclerosis, American Heart Association. Circulation 
1995;92(5):1355-74. 

[2] Rokitansky CV. A manual of pathological anatomy. Trans, by Day GE. Vol. 4. London: 

Sydenham Society, 1852. 
[3] Virchow R. Der ateromatose prozess der arterien. Wien Med. Wochenshr 1856;6:825- 

41. 

[4] Anitschkow N. Experimental arteriosclerosis in animals. In: : Cowdry EV, Editor. 

Arteriosclerosis - a survey of the problem. New York: Macmillan, 1933;27 1-322. 
[5] McGill HCJ, Geer JC, Strong JP. Natural history of human atherosclerotic lesions. In: 

Sandler M, Boune GH, editors. Atherosclerosis and its origin. New York: Academic 

Press, 1963;39-65 

[6] Fox JA, Hugh AE. Localization of atheroma: a theory based on boundary layer 

separation. Br Heart J 1966;28:388-99. 
[7] Fry DL. Acute vascular endothelial changes associated with increased blood velocity 

gradients. Circ Res 1 968;22 : 1 65-82. 
[8] Caro CG, Fitzgerald JM, Schroter RC. Arterial wall shear and distribution of early 

atheroma in man. Nature 1 969;223: 1 1 59-6 1 . 
[9] Haust MD, More RH. Development of modem theories on the pathogenesis of 

atherosclerosis. In: Wissler RW, Geer JC, editors. The pathogenesis of atherosclerosis. 

Baltimore: Williams and Wilkins, 1972; 1-1 9. 



Analytical Methods of Atherosclerosis Research 



309 



[10] Ross R, Glomset J. Atherosclerosis and arterial smooth muscle cell. Science 1973; 
180:1332-9. 

[11] Benditt EP, Benditt JM. Evidence for a monoclonal origin of human atherosclerotic 

plaques. Proc Nat Acad Sci USA 1973;70:1753-6. 
[12] Stemerman MB. Haemostasis, thrombosis and atherogenesis. Atherosclerosis Rev 

1979;6:105-46. 

[13] Texon M. Hemodynamic basis of atherosclerosis. Washington: Hemisphere, 1980. 
[14] Friedman MH, et al. Arterial geometry affects hemodynamics: a potential risk factor for 

atherosclerosis. Atherosclerosis 1983;46:225-31. 
[15] Solberg LA, Strong JP. Risk factors and atherosclerotic lesions: a review of autopsy 

studies. Arteriosclerosis 1983;3:187-98. 
[16] Debakey ME, Lawrie GM, Glaeser DH. Patterns of atherosclerosis and their surgical 

significance. Ann Surge 1985;201:115-31. 
[17] Giddens DP, Bharadvaj BK, Ku DN. Cerebral and peripheral hemodynamics. In: 

Altobelli SA, Voyles WF editors. Cardiovascular ultrasonic flowmetry. New York: 

Elsevier, 1985;125-145. 
[18] Schwartz CJ, et al. The pathogenesis of atherosclerosis: an overview. Clin Cardiol 

1991;14:1-16. 

[19] McGill HC Jr. Atherosclerosis: a moving target-key note presentation. In: Wissler RW 
editor. Atherosclerotic plaques advances in imaging for sequential quantitative 
evaluation. NATO SAI Series. New York: Plenum Press, 1991;219:1-9. 

[20] Grundy SC. Role of low-density lipoproteins in development of coronary artery 
atherosclerosis. In: Kreisberg RA, Segrest JP, Editors. Plasma lipoproteins and 
coronary artery disease. Cambridge: Blackwell Scientific, 1992:93-124. 

[21] Nerem RM. Vascular fluid mechanics, the arterial wall, and atherosclerosis. AS ME J 
Biomech Eng 1 992; 1 14:274-82. 

[22] Harry L. Goldsmith H.L and Takeshi Karino. Flow and vascular geometry. In: Hwang 
Ned HC, Turitto VT, Yen M RT. Editors. Advances in Cardiovascular Engineering. 
NATO ASI Series. New York: Plenum Press, 1992;235:127-51. 

[23] Fuster V. Human lesion studies. Ann N Y Acad Sci 1 995;748:207-24. 

[24] Berliner JA, Navab M, Fogelma AM. Atherosclerosis: basic mechanisms. Circulation 
1995;91:2488-96. 

[25] Breslow JL. Mouse models of atherosclerosis. Science 1 996;272:685-8. 

[26] Davies MJ. Stability and instability: two faces of Coronary atherosclerosis. The Paul 

Dudley White Lecture. Circulation 1996;94(8):2013-20. 
[27] Murray CJ, Lopez AD. Global mortality, disability, and the contribution of risk factors: 

Global burden of disease study. Lancet 1997;349: 1436-1442. 
[28] Hansson GK. Cell-mediated immunity in atherosclerosis. Curr Opin Lipidol 

1997;8:301-11. 

[29] Ciofalo M, Collins MW, Hennessy TR. Nanoscale fluid dynamics in physiological 

processes: A review study. Boston: WITPtess, 1999. 
[30] Sun Z, Larson IA, Ordovas JM, Barnard JR, Schaefer EJ. Effects of age,' gender, and 

lifestyle factors on plasma apolipoprotein A-TV concentrations. Atherosclerosis 

2000;151:381-388. 

[3 1 ] Faergeman, O. The revised joint guidelines. Atherosclerosis Suppl 2000; 1 :3-7. 
[32] Lusis AJ. Atherosclerosis. Nature 2000;407:233-41 . 



3 1 0 Harris H. Wang and Xing F. Wang 



[33] Risch NJ. Searching for genetic determinants in the new millennium. Nature 
2000;405:847-56. 

[34] Fayad ZA, Fuster V. Characterization of atherosclerotic plaques by magnetic resonance 

imaging. Ann NY Acad Sci 2000;902:173-186. 
[35] Wang HH. Analytical models of atherosclerosis. Atherosclerosis 2001;159:1-7. 
[36] Kruth HS. Lipoprotein cholesterol and atherosclerosis. Curr Molecular Medicine 

2001;1:633-53. 

[37] Libby P. Inflammation in atherosclerosis. Nature 2002;420:868-74. 

[38] World Health Report. Report of the Director General, WHO. Geneva: 1998. 

[39] American Heart Association. 2002 Heart and stroke statistical update. Dallas: 

American Heart Association, 2001 . 
[40] Shephered J, Betteridge DJ, Durrington P. Strategies for reducing coronary heart 

disease and desirable limits for blood lipid concentrations: guidelines from the British 

Hyperlipidaemia Association. Br Med J 1987;295:1245-6. 
[41] Study group of the European Atherosclerosis Society. The recognition and management 

of hyperlipidaemia in adults. A policy statement of the European Atherosclerosis 

Society. Euro Heart J 1988;9:571-600. 
[42] Canadian lipoprotein conference at hoc committee on guidelines for 

dyslipoproteinemias. Guidelines for the detection of high risk lipoprotein profiles and 

the treatment of dyslipoproteinemias. Can Med Assoc J 1990;142:1371-82. 
[43] National cholesterol education program. Second report of the expert panel on detection, 

evaluation, and treatment of high blood cholesterol in adults (adult treatment panel II). 

Circulation 2002; 106(25): 1333-445. 
[44] National Center for Health Statistics. National health and nutritional examination. 

survey (III): 1994. 

[45] Ross, R. Atherosclerosis - an inflammatory disease. N Engl J Med 1999;340(2): 1 1 5-26. 
[46] Allard C van der Wal, Onno J. De Boer and Anton E. Becker. Inflammation in coronary 

atherosclerosis-pathological aspects. In: Mehta JL. editor. Inflammatory and infectious 

basis of atherosclerosis. Basel: Birkhauser Verlag, 2001 ;23-47. 
[47] Li AC, Glass CK. The macrophage foam cell as a target for therapeutic intervention. 

Nature Medicine 2002;8(1 1): 1235-42. 
[48] Binder CJ, et al. Innate and acquired immunity in atherogenesis. Nature Medicine 

2002;8(ll):1218-26. ' : 

[49] Beere PA, Glagov S, Zarins CK. Retarding effect of lowered heart rate on coronary 

atherosclerosis. Science 1984;226:180-2. 
[50] Kannel WB, Kannel C, Paffenbarger RSJ. Heart rate and cardiovascular mortality: The 

Framingham study. A m Heart J 1 987; 1 1 3 : 1 489-94. 
[51] Fry, D.L. In Cerebral-vascular disease: Hemodynamic forces in atherogenesis. New 

york: Raven Press, 1976;77-95. 
[52] Caro CG. Arterial fluid mechanics and atherogenesis. Clin Hemorheol 1 982;2: 131-6. 
[53] Ku DN, Giddens DP, Zarins CK, Glagov S. Pulsatile flow and atherosclerosis in the 

human carotid bifurcation - positive correlation between plaque location and low 

oscillating shear stress. Arteriosclerosis 1985;5:293-302. 
[54] Lei M, Kleinstreuer C, Truskey GA. Numerical investigation and prediction of 

atherogenic sites in branching arteries. J Biomech Eng 1995;17:350-7, 



Analytical Methods of Atherosclerosis Research 



311 



Hademenos GJ, Massiud TF. Biophysical mechanisms of stroke. Stroke 1997;28:2067- 
7 

Gimbrone MA Jr, Topper JN, Nagel T, Anderson KR, Garciacardena G. Endothelial 
dysfunction hemodynamic forces, and atherogenesis. Ann NY Acad Sci 2000;902:230- 
40. 

Fisher WR, Granade ME, Mauldin JL. Hydrodynamic studies of human low density 
lipoproteins. Evaluation, of the diffusion coefficient and the preferential hydration. 
Biochemistry 1971;10:1622-9 

Fung YC Biomechanics, Circulation, Second edition. New York: Springer, 1997. 
Adamson AW. Physical chemistry of surfaces. 6 th edition. New York: John Wiley, 
1997. 

Welty JR, Wicks CE, Wilson RE. Fundamentals of Momentum, heat and mass transfer. 
4 th edition. New York: John Wiley, 2001 . 

Nichols WW, O'Rourke MF. McDonald's Blood in Arteries. 3rd edition. Philadelphia: 
Butler & Tanner, 1990. 

Schlichting H. Boundary layer theory, 6 th edition. New York: McGraw-Hill, 1968. 
Gutstein WH, Schneck DJ, Marks JO. In vitro studies of local blood flow disturbance in 
a region of separation. J Atherosclerosis Res 1968;8:381-8. 

Davies PF, et al. Turbulent fluid shear stress induces vascular endothelial turnover in 
vitro. Proc Natl Acad Sci USA 1986;83:2114-17. 

Liepsch DW. Effect of blood flow parameters on flow patterns at arterial bifurcations- 
studies in model. In: Liepsch DW, editor. Blood flow in large arterials: applications to 
atherogenesis and clinical medicine. Munich: Karger, 1990;63-76. 
Walburn FJ, Sabbah HN, Stein PD. Flow visualization in a model of an atherosclerotic 
human abdominal aorta. J Biome Eng 1981;103:168-70. 

Motomiya M, Karino T. Flow patterns in the human carotid artery bifurcation. Stroke 
1984;15:50-6. 

Could AL, et al. Cholesterol reduction yields clinical benefit: Impact of statin trials. 
Circulation 1998;97:946-52. 

National cholesterol education program expert panel. Report of the national cholesterol 
education program expert panel on detection, evaluation and treatment of high blood 
cholesterol in adults. Arch Intern Med 1988;148:36-69. 

Bargeron CB, Hutchins GM, Moore GW. Distribution of the geometric parameters of 
human aortic bifurcations. Atherosclerosis 1986;6:109-13. 

Tunito VT, Baumgartner HR. Platelet deposition on subendothelium exposed to 
flowing blood: mathematical analysis of physical parameters. Trans Am Soc Artif 
Organs 1997;21:593-600. 

Diana JN. Tobacco smoking and atherosclerosis: overview. In: Diana JN, editor. 
Tobacco smoking and atherosclerosis. New York: Plenum Press, 1990;l-7. 




Atherosclerosis 159 (2001) 1-7 



ATHEROSCLEROSIS 



ELSEVIER 



www.elsevier.com/locate/atherosclerosis 



Review article 



Analytical models of atherosclerosis 



Harris H. Wang * 



' 1 



3W Consulting Company, 61 Houghton St. No. 5, Worcester, MA 01604, USA 
Received 8 December 2000; accepted 4 June 2001 



Abstract 



Atherosclerosis is responsible for ^50% of all mortalities in the USA, Europe, and Japan. An innovative approach for 
investigating atherosclerotic lesions using a thermodynamic model and a boundary value model and an analytical example of a 
human abdominal aorta vessel with an initial lesion are presented in the paper. Analytical results given by both models propose 
that increased transient boundary layer thickness and reduced surface energy of adhesion in regions of the arterial branch points 
play crucial roles in initial lesion formation. This study also improves the understanding of atherosclerotic mechanisms and helps 
to interpret clinical and experimental results. © 2001 Elsevier Science Ireland Ltd. All rights reserved. 
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Nomenclature 

L a axial adhesive length (cm) 

z h axial distance between point being studied and starting point of transient boundary layer (cm) 

axial shear stress of plasma fluid (g/cm s 2 ) 

/ axial transient inertial force (g/cm 2 s 2 ) 

f h axial transient inertial force per unit area at border of transient boundary layer (g/cm s 2 ) 

u z axial velocity of plasma fluid (cm/s) 

u h axial velocity of plasma fluid at border of transient boundary layer (cm/s) 

F e body force vector (g/cm 2 s 2 ) 

A contact area between lesion and endothelium (cm 2 ) 

p density of liquid (g/cm 3 ) 

v e eddy velocity (cm/s) 

a inner radius of tube (cm) 

F A interfacial shear resistance (g/cm s 2 ) 

E k kinetic energy of lesion (g cm 2 /s 2 ) 

a) mean transient frequency (Hz) 

u mean velocity of blood flow in abdominal aorta (cm/s) 

P pressure vector (g/cm s 2 ) 

F{z) shear force (g/cm s 2 ) 
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£ sl surface energy of adhesion (g/s 2 ) 

y x surface free energy of liquid (g/s 2 ) 

y s] surface free energy of liquid-solid interface (g/s 2 ) 

y s surface free energy of solid (g/s 2 ) 

h(z) thickness of transient boundary layer (cm) 

u velocity vector of fluid (cm/s) 

u { velocity vector of lesion (cm/s) 

ju viscosity of liquid (g/s cm) 

V volume of lesion (cm 3 ) 



1. Introduction 

Atherosclerosis is characterized by the thickening, 
hardening, and loss of elasticity of the inner arterial 
walls. These changes in the walls can lead to the 
obstruction of coronary bloodstreams. The earliest 
atherosclerotic lesions, called fatty streaks, are simple 
deposits of low-density lipoproteins (LDL) on the aor- 
tic endothelium. These streaks occur at a universal age 
of ten, but progress to more advanced stages at differ- 
ent rates in different people. In the second stage, the 
lesions are called fibrous plaques because they contain 
not only more smooth muscle cells, but also more 
collagens and elastic tissues. The final stage is the 
complex lesions that are of thrombus formation with 
deposits of fibrins and platelets. The American Heart 
Association (AHA) committee has defined a new clas- 
sification of the human atherosclerotic phases from the 
fatty streaks to the complex lesions [1]. 

Since atherosclerosis is a major cause of mortality 
and morbidity in the world, a huge volume of research 
has been focused on the disease [1-18]. These research 
have been concerned mainly with dieting, environment, 
genetics, and immunity. A number of studies also at- 
tempted to identify the components and characteristics 
of lesions at various stages. Ross [2] contributed his 
response-to-injury hypothesis to the understanding of 
cellular and molecular mechanisms of atherosclerosis. 
These investigations indicated that hemodynamics 
might be closely related to atherosclerosis. 

Atherogenesis is related to nanoscale fluid dynamics 
and macromolecular transport at the arterial endothe- 
lium [3]. The behavior of atherosclerosis is time-depen- 
dent and the dominant phenomena may occur in the 
range of 0-100 A from the endothelium. The location 
of atherosclerosis is also associated with flow separa- 
tion and turbulence. These factors cause the study of 
the effects of hemodynamics on atherosclerosis to be 
rather complex. Early works on these effects focused on 
endothelial transportation using excised tissue. Fry 
[19,20] initiated studies on the effect of increased shear 
stress on the endothelium and showed that endothelial 



cells remodel their shapes when the direction of flow is 
changed. Caro et ah [21,22] found that lesions occur in 
areas experiencing low and fluctuating wall shear stress. 
Recent reports discussed the progress in the studies of 
the effects of hemodynamics on atherosclerosis [23-29]. 

Fung [28] summarized these hemodynamic effects 
and pointed out that the most intriguing feature of 
atherosclerosis is its appearance only at certain arterial 
branch points. It is known that many factors influence 
atherogenesis including dieting, environment, diabetes, 
high blood pressure, genetics, and exercise. These fac- 
tors, however, affect not only the special branch loca- 
tions, but also the entire body. Then why will only a 
few feet, in miles of blood vessels in the human body, 
be destined to have atherosclerosis? Why do only LDL 
deposit on some arterial walls and not others? What is 
so unique about these locations? Throughout the his- 
tory of the study of atherosclerosis, many researchers 
have been focusing on the understanding of these fun- 
damental natures of atherosclerosis. This paper pro- 
poses some answers to these fundamental questions. 

An innovative approach for investigating atherogene- 
sis is developed using a thermodynamic model and a 
boundary value model. Conclusions given by this study 
propose two new factors in atherogenesis: increased 
transient boundary layer thickness and reduced surface 
energy of adhesion in regions of arterial branch points, 
which play crucial roles in atherogenesis. 



2. Thermodynamic model 

2.1. Thermodynamic model of an atherosclerotic lesion 

The arterial endothelium is always in direct contact 
with the blood flow and is borne by hemodynamic 
forces. Endothelial cells have synthetic and metabolic 
capabilities and act as selectively permeable barriers for 
the transportation of macromolecules. The interfacial 
forces acting on the endothelium can change endothe- 
lial functions and structures including increased molec- 
ular permeability, LDL accumulation, and endothelial 
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cell damage. In this model, an initial lesion adhered on 
the endothelium is investigated using thermodynamics 
and energy conservation principles. 

Let us consider the blood flow through a branched 
circular arterial tube of internal radius a under an axial 
transient inertial force /, where (r, 0, z) is a set of 
cylindrical polar coordinates with the z-axis coinciding 
with the axis of the tube. Taking a small fluid element 
near the inner tube wall at a branch point where r = a 
and z = 0, the kinetic energy of the element E k is 
5'P'Vuf, where p is the density of the fluids V is the 
volume of the element, and w, is the velocity vector of 
the element. If an initial lesion occupies the element, E k 
becomes the kinetic energy of the lesion. This lesion is 
characterized by volume V, axial adhesive length of the 
lesion on the endothelium X a5 contact area between the 
lesion and the endothelium A, surface free energy y sl , 
velocity vector u u and interfacial shear resistance F a 
that is the force resisting the movement of the lesion 
(Fig. 1). All symbols of the paper are denoted in the 
nomenclature. 

2.2. Energy equation 

The adherence of an initial lesion on a plasma-en- 
dothelial interface can be treated as a thermodynamic 
problem in a mixing system because thermodynamics 
deals with macroscopic and statistical behaviors of 
interfaces rather than with details of their molecular 
structures. Surface energy of adhesion E sb one of ther- 
modynamic terms, is the energy that is required to 
bring LDL from the interior of the plasma fluid onto 
the arterial endothelium and form an interface. Blood 
flow is to obey the principles of conservation of energy, 
mass, and momentum. 

Consider a lesion under an axial transient inertial 
force / (Fig. 1). Based on thermodynamics and energy 
conservation principles, the energy equations for the 
lesion are created as 

E* = \f^L<A % (1) 

E* = \f*L*> (2) 
E* = E*A 9 (3) 



where £ sl = y s + y x - y sl [30]. 

Eq. (1) describes the equivalence between the kinetic 
energy of a lesion E k and the work done by the total 
interfacial shear resistance F. d A along the adhesive 
length L a . Eq. (2) represents the equivalence between 
the surface energy of adhesion E sl and the work done 
by F a at a liquid -solid interface. Substitution of Eq. (2) 
into Eq. (1) yields Eq. (3). 

If E k is greater than E sl -A, then the lesion will not 
form because the kinetic energy of the lesion E k over- 
comes the total energy of adherence E si -A before the 
adherence of the lesion can occur. 



3. Dynamic boundary value model 

3.1. Governing equation 

Studies of hemodynamics [21-23,28,31-35] suggest 
that atherosclerotic lesions are located in regions of 
reduced shear stress of the fluid, which is often associ- 
ated with flow separation and turbulence. A transient 
boundary layer, led by the flow, occurs in these regions. 
The studies also indicated that blood flow near inner 
arterial walls consists mainly of plasma fluids. 

In the dynamic model, the plasma viscous flow in the 
transient boundary layer near the inner tube wall at an 
arterial branch point is investigated because lesions 
arise from this layer. Since plasma fluid is known to be 
Newtonian, the flow is governed by the Navier- Stokes 
equation. The following restrictions are imposed to 
simplify the dynamic model. 

L The plasma may be treated as an incompressible, 

homogenous fluid. 
2. The flow is laminar. 

Let us consider an incompressible laminar flow 
through a circular tube. The governing equation is 

P'«^ e -gradP + /rV 2 ' W , . (4) 

where u is the velocity vector of the fluid, F e is the body 
force vector, P is the pressure vector, p is the viscosity 
of the fluid, p is the density of the fluid, and /rV 2, w is 
the viscous force. 
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3.2. Assumptions 

The following assumptions are imposed to simplify 
the boundary value model. 

1 . The gravitational potential in the fluid has a negligi- 
ble effect because the transient boundary layer is 
thin. 

2. The rate of flow is steady, and there is a constant 
pressure flow. 

3. The flow does not slip at the inner wall of the tube. 

4. The length of the tube is long compared to the 
region being studied, so that edge effects become 
negligible. 

5. The diameter of the tube does not vary with internal 
pressure. 

6. There is no convection in the flow. 

These assumptions and previous restrictions are valid 
in most studies of hemodynamics [28,35,36]. 

3.3. Simplified governing equation 

Let us use the same notations and consider an incom- 
pressible plasma flow in a circular arterial tube with 
ranges {a - h{z)) ^r^a, 0^0 ^2n, and — oo < z < 
oo, where h(z) is the thickness of the transient 
boundary layer. Using cylindrical polar coordinates (r, 

0, z), we yield /rV 2 -« = ^*j^ r '~[~^ wnere "z is a 

function of r only. 

Following the above assumptions, Eq. (4) is sim- 
plified to 



and 



-ft- 



d« 2 
d7' 



(5) 



(6) 



where u z is the axial velocity component of the plasma 
fluid, T re is the axial viscous shear stress of the fluid, 
and / is the axial transient inertial force. 

3.4. Boundary condition 

According to previous assumptions, the boundary 
conditions of the model are 



u z = 0 at r = a, 
z TZ = F(z) at r = a, 



(7) 
(8) 



where F(z) is the shear force acting on the inner wall of 
the tube. 

If a fatty streak appears at r = a and z = 0, we can 
substitute F(z) — F a into Eq. (8) to yield 



where the interfacial shear resistance F a is defined in the 
previous model. Through Eq. (9), we can now join the 
thermodynamic model and the boundary value model. 

3.5. Compatibility condition 

As stated previously, there is often separation and 
turbulent flow at the location of atherosclerosis. The 
flow at these sites consists of two different types. The 
first is the plasma flow in the transient boundary layer 
that is governed mainly by the viscous force. The 
second is the blood flow in the rest of the regions of the 
tube that is governed by the axial transient inertial 
force. Hence, the viscous force and.the transient inertial 
force must equal at the border of the layer for a unique 
analytic solution. The compatibility condition is created 
as 

T rz =/ h atr = 0-/2(z), (10) 

where f h is the axial transient inertial force per unit 
area. We may take f h = f(a — h{z)). 

3.6. Analytical solutions 

Eqs. (5) and (6) combined with Eqs. (7), (8) and (10) 
yield the axial shear stress of the plasma fluid T ra , the 
axial velocity component u z , and the compatibility 
equations. (Detailed derivations are omitted due to 
limited space.) 



F(z)=f-[a-h(z) + 



2a 



Substitution of F(z) = F a into Eq. (11) yields 



Substitution of Eq. (2) into Eq. (14) yields 
*(z) 2N 



oo 

(12) 
(13) 

(14) 
(15) 



3. 7. The effects of turbulence eddies on atherosclerosis 

Turbulence eddies do not originate in the transient 
boundary layer, however, these eddies may enter the 
layer from the side r <{a — h(z)) [37,38], If we consider 
the effect of the eddies in the layer on atherosclerosis, 
the thickness of the layer h(z) in the region near the 
branch point can be expressed approximately as 



T re = F a at r = a, z = 0, 



(9) 



h(z) « B-v c 



(16) 
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Fig 2 The relationship between the axial velocity of the plasma fluid 
2 the thickness of the transient boundary layer. 
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4. Conclusions 



,^ tt- o 1 we find that the shear 
(l)FromEq • 6), Fig.2 and 3 _ wet ^ ^ 

stress of the P^ ^/l, and inversely propor- 
surface energy o adhesioj ^ d layer 

tional to the thickness ^j"™'" in increa sed 
h{r) This means that reduced x n results m 

h ^VlZi El (16) indicate that increased 

(2) t-ig. J anu ^4 v / d ed d y velocity v e . 
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t lacto T*e formation of initial lesions 
a domtnant facto, m m teduced 

because, based on botl , ol tn * P is more easi i y . 

A" a— asis occurs ,„ re^ou of reduced 
t„. According to Con. 1 and tq. V * 

pSU-lndo.he.ia, interface and iniua. lestona form 
more easily in regions of reduced 
(4) Clinical «-^^^J£lii w.tb 

Eq (16), lesion formation occurs at the * ioc 

r^r?Sd^3S- 

reduced interfacial shear resistance ^' a * d a 

endothelial defects initially atherogenesis is 

(6) Clinical results [3-5] sho* ^ g , 
caused by depositions of LDL on tne arte 



4.1. Analytical example 



transient frequency o> - £ Hz, densrty , S ^ 
plasma viscosity M - L2g ■ c > rf ^ 

t -0 2 cm, and the range 01 ui^ t 

e cm - 2 s " 2 . Substitutions of a = 0.5 cm, p ^ 8 
-i / ,n ofm - ! s- 2 r = 0.4 cm into Eq. (12) yield 
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Fig. 3. The relationship between 
the thickness of the transient boundary layer. 
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Hum. According to the thermodynamic model and Eq. 
(3), LDL rather than HDL (high-density lipoproteins) 
are drawn from the interior of the plasma fluid onto the 
plasma-endothelial interface because the lower density 
of LDL yields a lower E sl than the £ sl for HDL 
deposition. Therefore, reduced £ s , and less energy re- 
quired for LDL deposition leads to atherogenesis. 

(7) Clinical results [4,5,28,35] show that initial lesions 
do not appear in veins and capillary vessels even 
though there is a lower r rz in veins than there is in 
arteries. This phenomenon can be explained by the lack 
of a transient boundary layer in veins and capillary 
vessels. The presence of a transient boundary layer is 
important to initial lesion formation. 



5. Mechanism of atherosclerosis 

The above conclusions propose that locations of 
atherosclerotic lesions are associated with regions of 
increased transient boundary layer h(z) and reduced 
surface energy of adhesion E sb which lead to atherogen- 
esis in only a small region of the vast systems of blood 
vessels. Analytical results of both models predict that 
increased eddy velocity v e and reduced E sl in these 
regions of the arterial branch points are main causes for 
the formation of initial lesions. 

More specifically, increased h(z) in regions of the 
arterial branch point is responsible for increased v c and 
reduced viscous shear stress T ra . Reduced E sl and in- 
creased v e contribute to bring LDL from the interior of 
the plasma fluid to the plasma-endothelial interface, 
which lead to the formation of a lesion. The LDL then 
passes through the arterial endothelium and stimulates 
monocytes to enter this area. The monocytes may then 
engulf cholesterol, and the lesion progresses to a more 
advanced stage. 

6. Discussions and prospects 

In this paper, the viscous flow of the plasma fluid in 
an arterial vessel is characterized by using symbols w 2 , 
r^, h(z), u e , a,f, y b //, and p. The plasma-endothelial 
interface is characterized by using symbols E sh y s , y sb 
F a , L a , and A. An initial lesion is characterized by using 
symbols p, V, E k , E sb y sb u h A, and L a . The interfacial 
shear resistance is characterized by using symbols F z , 
L a , and A, The interactions among the plasma flow, 
interfacial shear resistance, and endothelial cells at a 
plasma-endothelial interface are investigated using 
Eqs. (1)— (16). Analytical results show that atheroscle- 
rosis and endothelial defects are closely dependent on 
h{z), E sb t? e , and F a . Increased thickness of the transient 
boundary layer h(z) and reduced surface energy of 
adhesion E sl are controlling factors in atherogenesis. 



The conclusions given by this study offer some in- 
sight to the understanding of the fundamental behav- 
iors of atherosclerosis. These conclusions are supported 
by clinical and experimental results [3-5,21-23,28,31- 
35], which indicate that the analytical method used is 
successful as a potential tool in the study of atheroscle- 
rosis and other human lesions. This approach may help 
studies in the effects of interfacial shear resistance F a 
and shear stress of the fluid on the molecular 
biology of endothelial cells, LDL accumulation, molec- 
ular permeability, stress-sensitive gene expression, and 
signal transduction. 

Future works also include: 

1. Investigating the intermediate lesions in transition 
between the fatty streaks and fibrous plaques using 
analytical methods because these lesions represent 
the critical stage in the acceleration of 
atherosclerosis. 

2. Analyzing the effects of turbulence eddies on 
atherosclerosis because the transfer of macro- 
molecules from the interior of the blood fluid into 
the arterial endothelium are closely dependent on 
these eddies. 

These works will enable us to study environmental, 
dieting, lifestyle, high blood pressure, exercise, and 
genetic factors associated with atherosclerosis when the 
reward of prevention and the potential for regression is 
maximal. 
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